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Recently there has been a surge in interest in using video-microscopy techniques to infer the local 
mechanical properties of disordered solids. One common approach is to minimize the difference 
between particle vibrational displacements in a local coarse-graining volume and the displacements 
that would result from a best-fit affine deformation. Effective moduli are then be inferred under 
the assumption that the components of this best-fit affine deformation tensor have a Boltzmann 
distribution. In this paper, we combine theoretical arguments with experimental and simulation 
data to demonstrate that the above does not reveal information about the true elastic moduli of 
jammed packings and colloidal glasses. 


I. INTRODUCTION 

Characterizing the elasticity of soft disordered materi¬ 
als is challenging, in part because these systems often lie 
at the boundary of where classical theories of elasticity 
are applicable In jammed harmonic sphere packings, 
for example, a length scale below which continuum elas¬ 
ticity breaks down has been explicitly identified, and this 
length diverges upon the approach to the jamming transi¬ 
tion [5] . Another complication derives from the presence 
of nonaSine distortions that disordered materials experi¬ 
ence in response to imposed deformations. It has been 
argued that these nonaffinities arise from spatial fluctu¬ 
ations of local elastic moduli [3] and finite-temperature 
effects [11 IS]. 

In the last decade there has been a surge of experi¬ 
ments and simulations that aim to calculate elastic con¬ 
stants from particle-level fluctuations [Mil]. One com¬ 
mon approach focuses on thermally induced microscopic 
strain fluctuations I1IH1II1II3I. In these studies, the fluc¬ 
tuations of a locally-defined strain field are aggregated 
over time to arrive at a distribution of strains at each 
point. It is assumed that these strains are drawn from 
a Boltzmann distribution whose weight defines a local 
elastic modulus. Global elastic properties of the system 
are then computed by averaging or coarse graining these 
local moduli. In the case of crystalline systems mis] 
it is relatively straightforward to construct a local strain 
field on a per-particle basis by appealing to an unde¬ 
formed lattice. In the case of amorphous systems, how¬ 
ever, such an identification is no longer possible. Instead, 
it has been suggested that one should construct a coarse 
grained strain field by computing the best affine approx¬ 
imation to the collective motion that some neighborhood 
of particles undergo [5|. 

In this paper we argue that thermally induced particle 
motion in amorphous solids does not permit the use of lo- 
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cal affine strain distributions (computed as best-fit affine 
transformations) to compute elastic moduli. To demon¬ 
strate this concept we consider two systems: a simula¬ 
tion of harmonic disks in two- and three-dimensions and 
a quasi-two-dimensional experimental colloidal system, 
studied and described in detail in Ref. |T4|. We rely, in 
particular, on simulations of harmonic disks since their 
elastic moduli have been very well characterized by inde¬ 
pendent measurements and theoretical analyses. More¬ 
over, these systems can be simulated at arbitrarily low 
temperatures to ensure that we are truly in the regime 
where linear response is valid. One quantihcation of the 
degree of structural correlations captured by the best- 
fit affine transformation. A, are their dependence on the 
coarse-graining scale L. If A is to relate to linear elastic¬ 
ity then we should expect var[Aa^] ^ L~'^ as expressed 
by HI, where a and /3 specify the indices of the deforma¬ 
tion tensor. We will use this simple scaling relation as 
a benchmark throughout our analysis, and we will show 
that the structural correlations captured by A^/j are sys¬ 
tematically too weak to relate to the elastic moduli. 

We wish to emphasize that defining best-fit local 
strains in actively deformed materials can, in conjunc¬ 
tion with knowledge of the local stresses, still lead to 
meaningful information about the elasticity of disordered 
materials m- Additionally, tracking particle positional 
information over a long time can be effectively used to 
estimate the covariance matrix of the system El- How¬ 
ever, we argue that if there are only thermally induced 
fluctuations, then the distribution of best-fit strains does 
not contain information about the system elasticity. 

In Sec. |ll|we introduce the formalism commonly used 
to extract affine strains from thermal fluctuations, and 
we highlight some problems with interpreting exponential 
fits of the associated distribution to obtain elastic moduli. 
In Sec. Ill we study the simulations in a regime where the 
correlations between particle positions and displacements 
are relatively small, and we show that a simple statisti¬ 
cal model completely describes the distributions of local 
strains. The statistical model accurately predicts the dis- 
triubtions for for all temperatures, coarse-graining sizes, 
and pressures, and we show that these measurements are 
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emphatically not connected to the elastic moduli of the 
systems. In Sec. IV we discuss simulation and experi¬ 
mental measurements in a regime with increasingly large 
particle position and displacement correlations. While 
the correlations that enter the calculation of the best- 
fit affine deformation tensor in this regime might be ex¬ 
pected to enable one to deduce the elastic moduli, we 
again find that this is not the case. We discuss these re¬ 
sults and their consequences for interpreting experimen¬ 
tal data in Sec. |Vj The Appendices present a reformula¬ 
tion of the quantity used to define Aq,^ (Appendix 

0 , the statistical model used to understand the data in 
Sec. |III| (Appendix!^, and the details of the simulations 
and experimental systems studied (Appendix [C|). 


II. IDENTIFYING LOCAL STRAINS FROM 
THERMAL FLUCTUATIONS 

We begin by introducing the formalism that has most 
commonly been used to extract affine strains from ther¬ 
mal fluctuations. As mentioned above, the most-studied 
practical solution to the problem of nonaffinities is to find 
the best-fit affine transformation of local particle posi¬ 
tions at time t — At onto particle positions at time t, 
and then study distributions associated with this affine 
component of particle motion. This scheme is commonly 
done using a measure of nonaffinity, originally in¬ 

troduced by Falk and Langer [16] . One version of this cal¬ 
culation considers a local square or cubic coarse graining 
volume of side length L centered at point R. The mo¬ 
tion of the particles j in that local volume are tracked 
between times t — At, and t and then one computes the 
deviation of their displacements from those described by 
a best-fit affine transformation over that time window 
[TB] . Explicitly, 


U2(t,At)=^^fr“(t)-i?“ 

j a ^ 

+ Aafs) r^j{t-At)-R^ 


( 1 ) 


Here the Greek indices run over the Cartesian coordi¬ 
nates, fj{t) is the position of particle j at time t, A is an 
affine transformation tensor, and Sap is the Kronecker 
delta. This quantity is then minimized over all possible 
affine transformation tensors, Aq,^: 

At) = minD^it, At). (2) 

Ac^i 


Solving for the minimizing affine transformation is 
straightforward. Defining 

Xc.p = - i?“) X (rf (t - At) - RP), (3) 
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Y^p = Y^irfit - At) - i?“) X (rf (t - At) - R^), (4) 
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the best-fit tensor can be written as 

Ko^p = ( 5 ) 

7 

The standard approach [5| has then been to assume that 
these strains are drawn from a Boltzmann distribution 
whose energy is given by the elastic energy Ejp, = A^^L'^ 
where p, is the elastic constant. If this assumption holds, 
then the probability distribution of the squared strain 
components is given by H(A^^) ~ exp(—^A^^L'^/fcT) 
and the modulus can be extracted by fitting the loga¬ 
rithm of the probability distribution to a straight line. 

To test this basic claim we plot, in Fig. [T| the distribu¬ 
tions of A^^ as measured for small At for both 2D and 3D 
harmonically repulsive disks, as well as for the colloidal 
system at longer At. We find that this distribution has 
a pronounced curvature on a log-linear scale, suggesting 
that a simple exponential decay is a poor description of 
the distribution. A similar curvature in this distribution 
was reported by Rahmani et ah, where it was interpreted 
to result from a heterogeneous distribution of local mod¬ 
uli [17j . However, as shown in the figure, we find that 
the simulation distributions are accurately described by 
a simple distribution coming from the square of a sin¬ 
gle Gaussian random variable. Indeed, Ganguly et al. 
showed that a perfect hexagonal lattice at low temper¬ 
atures has a Gaussian distribution of A^y [5|, implying 
that there, too, the distribution of A^y would take a 
form and not appear as a straight line on a log-linear plot. 
Given that this qualitative signature of a heterogeneous 
distribution of elastic moduli can be completely repro¬ 
duced in the context of a simple statistical model with 
a single underlying Gaussian distribution (as discussed 
in more detail below), our results for disordered solids 
reinforce the message that it is hazardous to fit exponen¬ 
tial decays to these distributions and then interpret the 
result as elastic moduli. 

We note that an alternate definition of the best-fit 
affine deformation tensor that is commonly used, and is 
due to Cundall [H], first subtracts rigid body displace¬ 
ments of the local clusters center of mass before com¬ 
puting the affine distortion. Since amorphous materials 
can have large low-energy fluctuations compared to crys¬ 
talline systems, these center of mass displacements could, 
in principle, be important in our analysis, particularly 
for relatively small coarse graining lengths. However, we 
have carefully checked that the conclusions in our paper 
are insensitive to this definitional choice. 

In the following we will systematically study the dis¬ 
tributions of Aap from our simulations as a function of 
initial system pressure p, the time window At, coarse- 
graining scale L, and temperature T. We will supplement 
this with experimental data for a 2D colloidal sample. By 
varying the At at which we compute the best-fit affine 
deformation tensor between the ballistic regime and the 





3 



FIG. 1. (Color online) Probability distribution of L'^A.^y for 
simulations in d = 2 (upper red circles) and d = 3 (lower blue 
circles) of harmonic repulsive spheres at T = 10~®. Curves 
are best-fit distributions. Inset: Probability distribution 
of the unsealed for the experimental colloidal system. 


cage regime we can systematically tune the amount of 
correlation between displacements of a particle and its 
local environment. In no regime do we find a connec¬ 
tion between measured distributions of affine strains and 
elastic moduli. 


III. STRAIN MEASUREMENTS FOR SMALL At 


In the ballistic regime of small At (where time is mea¬ 
sured in units of r, the Lennard-Jones-like time unit of 
our simulation) there are only relatively small correla¬ 
tions between the frame-to-frame displacements of parti¬ 
cles and their initial positions. As such, we expect that 
the measurement of locally coarse-grained Aap can be un¬ 
derstood as the result of measuring single-point particle 
fluctuations and then averaging over a locally amorphous 
environment. Although our experimental data is well 
out of this regime, it is easily probed in our simulations. 
Clearly, in this regime the system lacks the particle- 
displacement correlations necessary to be described as 
a solid. Nevertheless, our exploration of this regime en¬ 
ables us to unambiguously identify an existing problem 
in how strain distributions have been analyzed in amor¬ 
phous solids [El EH- It also enables us to set up a conve¬ 
nient metric for how the strain variances must vary with 
coarse graining scale in order to be interpreted as moduli. 

As we show in Appendixj^ the calculation of the best- 
fit affine deformation tensor in a local coarse-graining 


volume can be usefully formulated as 

^ 7^17 j > (®) 

j \ 7 / 

where the first sum is over all particles j in the local 
coarse-graining volume, Aja is the frame-to-frame dis¬ 
placement of particle j in the a direction, the are 
quantities related to the initial positions of all particles 
in the local volume, and is the 7 component of the 
position of particle j at time t — At. In the limit of small 
correlations between displacements and local structure, 
then, Aa /3 can be approximated as a sum of random vari¬ 
ables, where Aja is drawn from a Gaussian distribution 
whose width, a a, is set by the density and temperature of 
the system and which is uncorrelated with the structural 
parameters 

In Appendix we combine this idea with the simplest 
possible model of the spatial structural parameters 
treating the positions of particles within the local coarse- 
graining volume as being uniformly distributed with no 
excluded volume. By the central limit theorem this sim¬ 
ple statistical model predicts that Aap will have an ap¬ 
proximately Gaussian distribution. Given a measure¬ 
ment of cta and the particle number density. Appendix 
[B] provides a prediction for its variance. 

Despite the naivete of this model, we find that the 
distributions of A^js when measured with a small At 
are remarkably well described by sums of Gaussian ran¬ 
dom variables multiplied by uncorrelated, ideal-gas-like 
structural parameters. Figure [^demonstrates this, show¬ 
ing that the distributions of A^j^, and Aj;y collapse when 
scaled by the appropriate powers of temperature and 
coarse-graining scale predicted by the model. Addition¬ 
ally, as expected by the model, the distributions for every 
component of Aq,^ is nearly identical. Even more remark¬ 
ably, the simple statistical model predicts the variance of 
the observed Gaussian distributions to within 10%. 

Given that the variances are so well-described by 
a statistical model with no positional correlations and 
no correlations between particle positions and displace¬ 
ments, any effort to extract elastic moduli from this 
measurement is doomed to failure. As a simple demon¬ 
stration of this failure, we follow Ganguly et al. and 
interpret var [Aa;^; -|- Ayy] as the bulk compliance and 
var [Axy + Ayx\ as the shear compliance [5] (given that 
P(A^^) is so well-described by a x? distribution it makes 
little sense to fit a straight line to it). Explicitly, as a 
function of coarse-graining volume L the relationship be¬ 
tween the bulk and shear modulus and local strain fluc¬ 
tuations is given by 

var [A.. + A,,] = ^ {B{L) + G{L))-^ (7) 

var[A.,+A,,,] = ^(G(L))-\ (8) 

We then plot the moduli - the inverse compliances - in 
Fig. ^ as a function of pressure for our simulated sys- 
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FIG. 2. (Color online) Collapsed probability distributions of P{h.xx) (blue points) and P{Axy) (red points) as computed 
for At = 2t. (A) Scaling collapse of P{Aap) for 2D simulations with ^/T for coarse-graining scale L = 2, p = 10~^ and 
temperatures of T = 10“®, 2 x 10“®, 4 x 10“®, 6 x 10“®, 8 x 10“®, 10“^. Dashed line is the prediction from Appendix [ b| 
(B) Scaling collapse of P{Aap) for 2D simulations with for T = 10“®, p = 10“^, and L = 2, 4, 6, 8, 10, 12, 14. Dashed 
line is the prediction from Appendixj^for the L = 6 data set. (C) Scaling collapse of P{Aap) for 3D simulations with ® for 
T = 10"®, p = 10"^ and L = 2, 4, 6, 8, 10, 12, 14. 


terns. Notably, the ratio of the bulk to the shear modu¬ 
lus is constant, whereas it is known that these jammed 
packings have a ratio that scales with the square root of 
the pressure, G/B ^ ^/p [19]. Both the bulk and shear 
moduli scale with the true bulk modulus of the system 
(the dashed line in the figure); the overall scale of fluctu¬ 
ations, a A, itself tracks the scaling of the bulk modulus. 
Thus, this measurement could be used to extract the scal¬ 
ing of the bulk modulus with, e.g., pressure, but not its 
absolute magnitude - of course, this scaling can be more 
more easily extracted by simply measuring the scale of 
the fluctuations directly. 

The fact that in the small At regime the Aap are single¬ 
particle quantities whose variance scales as L~'^ in two 
dimensions and in three dimensions is a useful ref¬ 
erence. Recall that the energy proposed by Schall et al. 
scales as E/p = and moreover that the mean of 

Aaj) is zero. It follows that if E/p is to be well-defined 
in the limit of large coarse graining sizes the variance of 
Aap must scale as L~'^. Therefore in the regime of small 
At, where var [Aq,^] ~ the variance of A^p cannot 

yield a well defined elastic modulus. Even though Fig. 
already demonstrated this to be the case, the necessary 
condition that var [A^^] ^ in the range of L con¬ 
sidered will be useful to keep in mind when we continue 
to the cage regime where we can no longer rely on an 
analytic model. 


IV. MEASUREMENTS IN THE PLATEAU 
REGIME OF THE MEAN-SQUARED 
DISPLACEMENT 

The above section is a useful illustration that distribu¬ 
tions of A^p should not always be interpreted in terms 
of elastic moduli. However, it is perhaps not surprising 



FIG. 3. (Color online) Inverse variance of A^x (blue circles) 
and individual particle displacement magnitude (red squares) 
for 2D simulations with T = 10“® and L = 4 as a function of 
pressure, normalized by the value at p = 10““^. Dashed line 
is the predicted scaling of the bulk modulus with pressure. 
Inset: Ratio of the variance of the diagonal to off-diagonal 
component of Aap as a function of pressure. The correspond¬ 
ing ratio of moduli scales as G/B ~ ^/p (dashed line). 


that the relatively uncorrelated nature of particle posi¬ 
tions and frame-to-frame displacements in the ballistic 
or crossover regimes of At can be understood in terms of 
sums of Gaussian random variables and not in terms of 
elastic moduli. By increasing the At we can continuously 
tune the degree of correlation between the particle dis¬ 
placements and the local structure, and one hypothesis 
is that it is precisely these correlations that allow one to 
probe the elastic moduli of the system in question. 
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In Fig. 1^ we show how the distributions P{Aais) 
change as At is continuously increased from the ballis¬ 
tic regime to deep into the plateau regime. The most 
apparent change is that the variance of the diagonal and 
off-diagonal components of Aap begin to separate. How¬ 
ever, we also see that the distributions become increas¬ 
ingly non-Gaussian. Figure highlights the change in 
the tails of these distributions, showing that as At is 
increased the distribution of P{A\p) continuously shifts 
from being extremely well-characterized by a fit to 
one with an apparent power-law tail. We note that the 
experimental data in this plot is much deeper into the 
plateau regime of the mean-squared displacement than 
our simulation data. Nevertheless, the robust presence 
of a power-law tail in both the simulations and experi¬ 
ments further emphasizes the danger of fitting exponen¬ 
tial decays to different parts of distributions and 

of interpreting those fits as elastic moduli. 

However, given that there are clearly some increasing 
correlations being picked up by the P(Aq,/ 3 ) distributions, 
one may wonder if the variances of these distributions are 
related to the elastic compliances. In Fig. [^we plot in¬ 
verse variances as a function of At and p. Once again 
we see that the ratio of moduli estimated in this way 
is independent of pressure, in stark contrast to the true 
elastic constants of these systems. Thus, although there 
are additional correlations in these data sets, they are 
not straightforwardly connected to the appropriate inte¬ 
grals over the covariance matrix that would allow one to 
correctly extract elastic constants. 

Finally, as in the case of our measurements at short At, 
we can attempt to understand with a finite-size analysis 
whether A^p has the necessary correlations with local 
structure to hope to extract elastic moduli from it. To 
this end, we once again consider how the variance of the 
strain tensor scales with the size of the coarse graining 
region used to construct it. In Fig. we see that the 
variance of Axx + Ayy (which one hopes to interpret as 
the bulk compliance) scales with L~^ in two dimensions, 
again in contrast to the L~'^ scaling that any sensible 
definition of strain must have. We conclude that al¬ 
though Aap has more correlation with structure in the 
cage regime than in the ballistic regime, it still cannot 
be used to give a well-defined elastic moduli. This scal¬ 
ing with L highlights a danger of fitting limited ranges 
of data without systematically checking the dependence 
on the coarse graining scale. In Appendix C, we explore 
the unphysical behavior of the inferred elastic moduli if 
one attempts to artificially fit small portions of the strain 
distributions to exponential decays. 


V. DISCUSSION 

With a combination of simulation and experimental 
data, we have demonstrated that for model disordered 
solids it is incorrect to connect exponential fits of P{A^p) 
with elastic moduli. In the limit of small At, we have pre- 
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FIG. 4. (Color online) (A) P{Axx) (red-to-orange color 
scale) and P{Axy) (dark-to-light blue color scale) for 2D 
simulation data with p = 10“^, T = 10“®, L = 4 for 
At/r = 2, 10, 20, 50, 100, 200. Data sets with smaller vari¬ 
ance correspond to shorter At. (B) P(A^y) normalized by 
the variance for the above At (red points correspond to the 
shortest At, light blue points to the longest At). Dark blue 
open circles are experimental data. The solid curve shows a 
unit-variance function, which matches the short-time data 
very well. 


sented an analytic model for the distributions of squared 
strain deformation tensor components. The model makes 
clear that the curvature of P{A\p) on a log-linear plot 
can be completely explained as a distribution coming 
from the square of a single underlying Gaussian distribu¬ 
tion, i.e., rather than from heterogeneous distribution of 
local elastic moduli. 

Furthermore, we have shown that for any choice of At 
the ratio of variances of A^x and Axy for these systems 
is essentially constant, with the ratio depending on the 
At window chosen but independent of the pressure of the 
sphere packings. In contrast, the global measurement of 
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FIG. 5. (Color online) (A) Inverse variance of P{Axx) as a 
function of At for pressures between p = 10“^ (blue, upper 
curve) and p = 10“^ (red, lower curve). Inset: Ratio of the 
variance of the diagonal to the off-diagonal components of 
Aq/s as a function of At (B) Ratio of variance of diagonal and 
off-diagonal components of Aq /3 as a function of pressure for 
At/r = 2,20,200. The corresponding ratio of moduli scales 
as G/B ~ 


G/B for these systems scales with the pressure of the 
packings, G/B ^ yPp. Thus, the ratio of the variances of 
these distributions does not correctly capture G/B, i.e., 
contrary to assumptions made in the literature [Hiin]. 
Moreover, we find that the variance of the strain in the 
cage regime scales with L~^ in two dimensions in both 
simulation and experiment. This implies that the strain, 
as computed via the best-fit affine transformation, does 
not have enough information about local structure to 
reliably extract elastic moduli. We contrast this with 
two-dimensional crystalline systems Pi where a simi¬ 
lar protocol reported variances that scaled with L~^ and 
elastic moduli consistent with other measurements were 



FIG. 6. (Golor online) Ya,r[Kxx. -\- Ayy], as a function of coarse 
graining size L in the simulation of two-dimensional harmonic 
disks for pressures from p = 10“^ (blue, bottom lines) to 
p — (red, top lines). Overlaid in black is a line with 

slope —3. Here Aap was calculated with At = 5 x 10^, deep 
into the caged regime. Inset: var[Ai,j; -|- Ayy\ as a function of 
coarse-graining scale (measured in microns) for area fractions 
<j} = 0.8625, 0.8695, 0.8822 (top to bottom). 


successfully extracted. 

Despite the fact that we find the variances of the best- 
fit affine strain distributions to be unrelated to the elastic 
moduli, we note that previous studies in amorphous ma¬ 
terials have found spatial strain-strain correlations with 
a quadrupolar signature in the nonlinear response involv¬ 
ing particle rearrangements [20j . This hallmark of con¬ 
tinuum elasticity is characterized by a power-law decay 
in the strain correlations far from the rearrangement. In 
the linear response of quiescent systems, by contrast, we 
are not aware of any experiment or simulation of amor¬ 
phous solids that shows that the best-fit affine strains 
arising from thermal motion have detectable power-law 
correlations, although such correlations must exist. In¬ 
deed, Rahmani et al. find that, in the absence of external 
strain, the strain-strain correlation function decays ex¬ 
ponentially (with a length scale on the order of a single 
particle diameter) [17]. 

The difference in this respect between crystalline and 
amorphous materials can be understood by considering 
the low-frequency excitations of the respective systems. 
In colloidal crystals (e.g., the hexagonal lattices studied 
in Refs. PP) the only vibrational modes are longitudi¬ 
nal and transverse sound modes that are (a) intimately 
related to the elasticity of the system and (b) spatially 
extended. We expect, via the equipartition theorem, that 
thermal fluctuations will populate these modes, lead¬ 
ing to extended and strongly correlated strain profiles. 
By contrast, disordered solids have a large population 
of vibrational modes that are extended but disordered 
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(e.g., in the boson peak). These modes span the sys¬ 
tem but have exponentially decaying local spatial cor¬ 
relations. Thermal population of these modes leads to 
strain profiles whose correlations are similarly exponen¬ 
tially localized; this is another argument for why mea¬ 
suring affine best-fit strains generated entirely by ther¬ 
mal fluctuations does not lead to meaningful information 
about the elasticity of disordered materials. In contrast, 
measuring these quantities in actively stressed or strained 
systems may still lead to meaningful information, as was 
discussed in a different context in Ref. [15] (see also be¬ 
low). 

The disordered solids whose elasticity we are testing 
have explicit length scales, whose scaling for our model 
system goes as I* ~ and ^ |5J [5T], be¬ 

low which continuum elasticity fails to describe the re¬ 
sponse of these solids to imposed forces and deforma¬ 
tions. We can estimate the magnitude of these length 
scales as ranging from I* « 7.1cr and Lt ~ 2.7a at 
p = 10“^ to I* « 71cr and Lt ~ 8.4(T at p = 10“"*^. The 
range of coarse-graining length scales considered may be 
compared to the characteristic size of structural hetero¬ 
geneities that lead to fluctuations in the local elastic mod¬ 
uli. In a jammed solid, these heterogeneities are expected 
to be on the scale of £t, or at worst £*. At the pressures 
studied, we have chosen coarse-graining scales that sat¬ 
isfy both L <^ £* and L ^ £*. In both of these limits, 
comparing these lengths with the data in Fig. [^suggests 
that the measurement of local affine strains is unable to 
detect the presence of these length scales, further under¬ 
mining the claim that var [A^p] is, on its own, intimately 
related to the local moduli. 

We note that other methods for defining local strain 
fields and then connecting these to local elastic moduli 
have been proposed. Tsamados et al. m studied a 
linear strain tensor, constructed from a continuous 
displacement field that was itself a coarse graining of lo¬ 
cal particle motion. In contrast to the methodology dis¬ 
cussed here, Tsamados et al. also measured a local stress 
tensor, computed by a similar coarse graining, and de¬ 
fined the local moduli to be the constant relating these 
coarse-grained strains to coarse-grained stresses. Unfor¬ 
tunately for experimental measurements, the computa¬ 
tion of the local stress fields requires knowledge of inter¬ 
particle forces, which are typically difficult to identify in 
experimental systems. However, it might be interesting 
to see whether the thermal fluctuations of this, or other, 
definitions of local strain can be used in an argument in 
the spirit of Schall et al.[5|. 

Another approach that has proven fruitful is to extract 
the particle-displacement covariance matrix from micro¬ 
scopic measurements. The bulk and shear moduli of the 
system can then be estimated from the inferred longi¬ 
tudinal and transverse speeds of sound in the material 
[14j . However, this method suffers from a few noted dis¬ 
advantages. Most significantly, a large amount of data 
is needed before the covariance matrix converges; this 
amount increases linearly with the number of particles in 


the system, L‘^, where L is the system length and d the 
dimensionality |22j . Second, disordered systems contain 
excess vibrational modes at low frequency that obscure 
the longitudinal and transverse acoustic branches of the 
phonon spectrum in systems that are too small; this effect 
scales as 1/L. Together, these issues limit the utility of 
the covariance-matrix-approach to systems that are nei¬ 
ther too large nor too small. This, then, explicitly limits 
the use of this tool when the distributions and spatial 
organization of elastic moduli are of interest. 

In closing, our results highlight the subtlety of measur¬ 
ing the elasticity in soft disordered systems: a method¬ 
ology that has been well-validated for two-dimensional 
crystalline systems fails spectacularly when applied to 
numerical simulations of disordered soft repulsive disks 
and laboratory experiments on colloidal packings. In 
light of this failure, we emphasize the critical importance 
of validating new methods of probing elastic constants 
by first testing them systematically against model sys¬ 
tems whose properties are known by more conventional 
elasticity measurements. 
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Appendix A: Defining local non-affinity 

In this appendix and the following we derive a simple 
statistical model that, for small At, almost completely 
captures the behavior of the components of the best-fit 
affine deformation tensors discussed above. For conve¬ 
nience, we first review an equivalent formulation of the 
^min language. In Appendix]^ we will employ this lan¬ 
guage to make simple estimates of the variances of Aap. 

1. Operator expressions for non-affinity 

Here we closely follow the language of Ganguly et al. 
[5]. In what follows roman indices will refer to parti¬ 
cles and Greek indices to spatial coordinates. We begin 
by defining the initial position of particle i, This 

could be the position of the particle at time t — At 
as in the definition, or we could take it to be 

the inherent structure position of particle i or its time- 
averaged position. Displacements from these initial po¬ 
sitions will be written as by Ui^{t) = ri^it) — For 
computing the local non-affinity for particle i we addi¬ 
tionally define displacements relative to that particle as 


Aj (t) = Uj {t) — Ui (t). Note that it is common to choose 
a reference position about which to define a local coarse 
graining volume, instead of a reference particle. In that 
case simply sets the origin of the local coordinate sys¬ 
tem, which does not change between time t— At and time 
t. In the following we will drop the explicit dependence 
on t in our expressions when it is clear from context. 

We next define an intensive measure of the local non¬ 
affinity of displacements relative to a reference particle 
in an analogous way to 




11(4 - A(r^ 

{ij) 


J)f 


(Al) 


where the sum runs over all particles j in the neighbor¬ 
hood considered, and N is the number of particles in that 
neighborhood. Without the factor of 1/N and taking 
to be the particle position at time t — At, the nonaffinity 
X is exactly equal to the definition of in the main 

text. Independent of the presence of the prefactor 1/N 
the tensor here is identical to the best-fit affine de¬ 
formation tensor defined in the introduction. Dividing 
out by the number of neighbors has been previously used 
to study thin films and pillars, where particles near the 
interface have many fewer neighbors than those in the 
center of the sample [ISIIH]. Nevertheless, since below 
we will be exclusively interested in the distribution of the 
components Aap the choice of an extensive or intensive 
definition of the total nonaffinity is irrelevant. 

In order to express both y, and especially A, in a con¬ 
venient operator form we define the following matrices. 
First, where d is the spatial dimension, we write the 
(lx dN) matrix 


Taking dx/dX and solving gives the minimizing affine de¬ 
formation: 

A = (R^R) R^A = QA. (A6) 


Given this minimizing A, the non-affinity can be written 
as 


X = ^ (A - RQAf 


= ^A'^ \1-2R{R'^R) ^R'^ + R{R^R) ' R^ 


, -1 


(A7) 

A 


= —A^PA. 
N 


The above expression defines a projection operator P = 
1 — RQ which projects components of A onto the space 
of non-afhne deformations. 


2. Specialization to two dimensions 


For concreteness, we explicitly write down an expres¬ 
sion for the components of A for a two-dimensional sys¬ 
tem. Taking a reference position to set the origin of 
our local coordinate system, we have Rja,jY = Sa-yrj^, 
for each particle j in the coarse-graining area that we 
choose. The matrix then has a simple struc¬ 

ture: 


where 




C B 0 0 \ 

B A 0 0 I 

0 0 C B ’ 

0 0 B A / 


(A8) 


A = (All, • ■ •, Aid, A 21 ,..., A2d, ■ ■ ■, Affd), (A2) 

which compactly writes all of the relative displacements 
in a convenient order. Next we define the (d^ x 1) matrix 

A = (All, • ■ •) Aid, • ■ • J Add)^ ■ (AS) 




A = 


ac—b^ ’ 


B = 


— b . 
ac—b’^ ’ 


(A9) 


C = 


ac—b"^ ■ 


Thus, in two dimensions the operator Q can be written 
as a combination of x d blocks, each of which looks 
like 


This simply unwraps the components of the best-fit A 
afhne-deformation tensor into a 1 dimensional array. Fi¬ 
nally, we define the {dN x d^) matrix 

~ (A4) 

This is a particularly convenient matrix with which to 
describe the initial relative positions of particles in the 
neighborhood of the reference particle. 

We are now in a position to express the non-affinity in 
a very compact fashion. With the above definitions we 
have 

X = [A - RXf (A5) 

= ^minA [A^A - A^RX - X^R^A + X^R^RX] 


Q = {R^Ry^R/^ 


I 

V 


Cr% + Brly 0 

Br% + Ar]y 0 

0 Cr° + By 

0 


\ 

/ 


(AlO) 

This lets us compactly write any component of A using 
A = QA, e.g. 


Axy — Ajx + Ar^y) ■ (All) 

j 

By writing the best-fit affine transformation tensor as a 
linear operator acting on the fluctuations it is already 
clear that one would not in general expect, e.g., an ex¬ 
ponential distribution of A^y for short At. The rela¬ 
tive displacements Aj can be assumed to be normally 
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distributed, after which the algebra of random variables 
suggests that has a ^ form. In the next section 
we show that a simple statistical model reproduces the 
distributions of Kxy that we observe in our simulations. 


Appendix B: Statistical model 


Here we show that in a disordered material we can 
use the algebra of random variables to accurately predict 
the distributions associated with at short At. As 
seen in Fig. and as could be anticipated from the 
functional form of Eq. |A11| in the absence of symmetry 
constraints and correlations, all of the components of A 
have nearly identical distributions when averaged over 
our disordered systems. Our goal in this section will be 
to predict the variance of as a function of the typical 
scale of the fluctuations of A and the size of the coarse- 
graining volume. To do so, we start from a simple model 
for single-particle positional distributions and build up to 
the distribution of the best-fit affine deformation tensor. 

For simplicity we focus on the two-dimensional case, 
and our dominant assumption will be a lack of structural 
order in the square coarse-graining cells. Hence, for a 
square coarse-graining square of side length L — 2R we 
take the to be uniformly distributed in {—R, R), i.e. 
to have a probability distribution given by 

F”.w={f |;;;|>^ ™ 


The building blocks of the best-fit affine deformation ten¬ 
sor involve sums of products of these single-particle dis¬ 
tributions. It is straightforward to show that 


Pj .0 ^0 (x) 

jx jy 



|x| < R^ 
|a;| > R? 


(B2) 


P(^^0J2{X) 


0 < a; < 

Ry/X 

0 otherwise 


(B3) 


To make further progress we invoke the central limit 
theorem to describe the a, &, and c random variables. Let 
n = p(2i?)‘^, where p is the number density, denote the 
average number of particles in a local coarse-graining vol¬ 
ume, and the symbol N (p., a) denote a Gaussian distri¬ 
bution with mean p, and width a. Then we approximate 


a{x) = Pfj .0 Y (x) « Af 


b(x) = P^q ^q (x) 
j=i 


'■N 


nR^ 2R^^Jn\ 



c{x), 


(B4) 


We next approximate the denominators that appear in 
the random variables A, B, and C, i.e. (ac — 6^). The 
part is trivial, and is given by 


Pb2 (x) 


3exp(^) 

\/2'KnR'^x 


(B5) 


The product ac can be written as the sum of two general 
distributions: 


ac = 


(a -I- cY 
4 


(q - cY 

4 


(B6) 



Note that both the second term in Eq. |B6| and the dis¬ 
tribution of the b^ have their weight centered about zero, 
whereas the first term in Eq. |B6| has a large positive 
mean. For simplicity we thus approximate the expres¬ 
sion (ac—6^) by a single non-central random variable: 


ac — b^ 



(B7) 


Defining aac = R^y/8n/45, the quantity 4(ac — b^')l<j\c 
is a non-central yY random variable with non-centrality 
parameter A = 5n/2 and number of summed normal vari¬ 
ables k = \. This allows us to write the first moment and 
variance of {ac — as 


(ac — Y) 


nR'^{2 + bn) 

' 45^"^ 


(B8) 


var [ac — 


8n2(l-h5n)i?® 

2025 


(B9) 


We now approximate the mean and variance of, e.g., 
A and B by the lowest order terms in the Taylor ex¬ 
pansion for the ratio of random variables, neglecting any 
covariance. That is, for random variables X and Y we 
approximate 



We find 


(41) 


var [A] 

{B) 

var[B] 


{X){YY , (A)var[r] , 

(BIO) 

(y)3 ' 

(y)3 ■ 

var[A](y)2 

{XY var[F] 

• (Bll) 

(U)4 

(y)4 + 

15 

120(1 -b5n) 

(B12) 

{2 + 5n)R^ 

(2-b5n)3i?2 

180 

1800(1 -f5n) 

(B13) 

n{2 + 5n)2i?4 

(2-t5n)4i?4 

0 


(B14) 

225 


(B15) 

n{2 + bnYR'^ 



The penultimate step is to consider the variance of the 
products AjxAxjy and Aj^Bcj^. We again simply as¬ 
sume that the are uniformly distributed and that the 
Aja are normally distributed with zero mean and width 
a A- Using the relation that the variance of a product of 
random variables is 


var [Ai • A 2 • • • A„] = [Q (var [A,] + (X^f) - n(^*)" 

i i 

(B16) 
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we have that 

var [Aj^Ar^y] 


var [Aj^Br°^] 

We are finally in position to evaluate the variance of 
the components of the best-fit affine deformation tensor. 
Since 


75cr 


A 






(2-h5n)2i?2 n(2-h5n)2i?2 

1800cri(l-h5n) 


(2 -t 5n)4i?2 


75ai 


n(2 -h 5 n) 2 i ?2 




(B17) 


(B18) 


Ajx {Btj^ + Arjy'j (B19) 

3 

we assume that the variance from each particle in the 
local coarse-graining volume contributes identically, and 
thus have 


var[Aa;y] « n (var [Ajj:Ar°y] + var [Aj^Br°^]) , (B20) 


where the variances in this equation are given by Eq. 


B17 A comparison between this equation and the simu¬ 
lation data (using the measured a a) is shown in Fig. 
where it is seen to be an excellent estimate of the vari¬ 
ance: with no adjustable parameters, and completely ig¬ 
noring correlations from excluded volume (or any other 
source), this simple model describes the variance mea¬ 
sured in the simulations to within 10%. 


Appendix C: System preparation 
1. Simulation details 

Our simulations are of frictionless packings with peri¬ 
odic boundary conditions composed of equal numbers of 
small and large spheres with a diameter ratio 1:1.4 and 
of equal mass, m. The particles interact with a repulsive, 
finite-ranged potential 



FIG. 7. (Color online) Mean-squared displacement (in units 
of particle diameters) for p — 10“^ '^, 10~^'°, 10“^ '*, 10“^'° 
(top to bottom). Inset. Mean-squared displacement 
in units of nm for the colloidal system for (j) = 

0.8625,0.865,0.8695,0.8775,0.8822 (top to bottom). 


the LAMMPS package [26]. For ease of comparing our 
timescales with the typical glassy crossover from ballis¬ 
tic to caged to diffusive behavior, mean-squared displace¬ 
ment curves for a subset of our 2D simulations are shown 
in Fig. 0 

In both the simulations and the experiments we com¬ 
pute the local non-affinity and best-fit affine deformation 
tensors as described in the text by partitioning the sys¬ 
tem into squares (cubes in 3D) of a given side-length. 
Thus, we compute with respect to a local origin of a 
coordinate system, rather than with respect to tagged 
reference particles. We have confirmed that this choice 
does not affect our conclusions. 


2. Experimental details 


= -o) (Cl) 

\ 0 r,j > aij, 

where is the distance between particles i and j, aij 
is the sum of the particles’ radii, and e determines the 
strength of the interaction. We report energies in units 
of e and distances in units of the average particle di¬ 
ameter. Time is measured in units of y/ejjma^. We 
used this model to study 1024-particle systems in 2D and 
4096 in 3D, for a range of pressures between p = 10“^ 
and p = 10“^. The initial configurations of these sys¬ 
tems were set by first placing the particles at random in 
an infinite-temperature configuration, and then quench¬ 
ing to T = 0 using a combination of linesearch meth¬ 
ods, Newton’s method, and the FIRE algorithm [^. We 
then perform low-temperature molecular dynamics using 


Our experimental systems are quasi-two-dimensional 
packings of poly(N-isopropyl acrylamide) (PNIPAM) mi¬ 
crogel particles. The full details of the experimental 
setup and data acquisition is reported in Ref. |17|. In 
brief, the disordered packing was prepared using a binary 
particle suspension with PNIPAM particles of two diam¬ 
eters: ai « 1.0 pm and a 2 ~ 1.4 pm. The sample was 
confined between two cover slips (Fischer Scientific) and 
then sealed from the edges with optical glue (Norland 
63) |2Hj . Since PNIPAM is a temperature-sensitive poly¬ 
mer, the particle diameters can be controlled by chang¬ 
ing the temperature. Thus, we tuned the effective pack¬ 
ing fraction of the sample, 4>, in situ using an objective 
heater (BiOptics). The temperature was set to a narrow 
range of 26.4 — 27.2°C so that the packing was above 
the jamming point. For each temperature studied, the 
trajectories oi N k, 4500 particles in the field of view 
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FIG. 8. (Color online) Normalized probability distributions 
of (Ajja; + Ayy)~2 for tho experimental colloidal system at a 
fixed volume fraction, truncated to a fixed dynamical range, 
for several different choices of the box coarse graining size. 
The fits are exponential decays, and the curves are calculated 
distributions for L = 1,3,4,11,20 microns. Inset. Inferred 
modulus from the exponential decay fits normalized by the 
measured value by the methods in Ref. [14] as function of L. 
Data sets are for (j) = 0.8625,0.865,0.8695,0.8775,0.8822. 


were extracted from a total of 30,000 frames of video us¬ 
ing standard centroid-finding and particle-tracking tech¬ 
niques mj. calculations were done by first sub¬ 

tracting the global drift of the sample and then using 
At = 0.273 s, a value which is well into the plateau re¬ 
gion of the mean-squred displacement [14j , as seen in the 
inset to Fig. 


Under experimental conditions, where the full distri¬ 
bution and its heavy tail may not always be accessible, 
it may be tempting to try to fit to an exponential decay 
to these distributions over some limited dynamic range. 
Thus, in Fig. we briefly mention the results of ap¬ 
proximating the observed experimental distributions of 
(Axx + Ayy)~'^ by exponential decays using an artificially 
restricted dynamic observation range. In the main plot 
we show exponential fits to the squared strain distribu¬ 
tions as a function of coarse graining size, and there is 
little agreement between these fits. In the inset we plot 
the bulk modulus that would be inferred from such fits 
relative to the bulk modulus as measured in Ref. [14]. 
There is no systematic trend suggesting that, in the large 
L limit where the method is ostesibly most sensible, the 
inferred modulus is asymptotically approaching the true 
value. This behavior could be anticipated from the re¬ 
sults shown in Fig. 
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